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Abstract 



We examine the hadronic interactions of quarkonia, focusing on the decays 
tp(2S) — ► J/ip-Kir and T(25) — ► Y(lS)irir. The leading gluonic operators 
in the multipole expansion are matched onto the chiral lagrangian with the 
coefficients fit to available data, both at tree-level and loop-level in the chiral 
expansion. A comparison is made with naive expectations loosely based on 
the large- iV c limit of QCD in an effort to determine the reliability of this 
limit for other observables, such as the binding of J/tj) to nuclei. Crossing 
symmetry is used to estimate the cross-section for inelastic irJ/ip — * mp(2S) 
scattering, potentially relevant for heavy ion collisions. The radiative decays 
ip(2S) — > J/-07r + 7r _ 7 and T(25) — ► Y(lS)7r + 7r _ 7 are determined at tree-level 
in the chiral lagrangian. Measurement of such decays will provide a test of the 
multipole and chiral expansions. We briefly discuss decays from the T(3S) 
and also the contribution from 7r's to the electromagnetic polarizability of 
quarkonia. 
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I. INTRODUCTION 



The strong interactions of quarkonia, such as the J / ip and T are much simpler than those 
of other hadrons as they are compact on the hadronic scale. Gluons with a wavelength A 
much larger than the "radius" , Tqq, of quarkonia have interactions suppressed in a multipole 
expansion by powers of Tqq/X. The leading operators involving two gluons correspond to 
interactions with external chromo-electric and chromo-magnetic fields with higher 

order interactions naively suppressed by additional factors of the quarkonium radius. In 
order to determine the interactions of quarkonia with other hadrons nonperturbative matrix 
elements of the gluonic operators between the hadronic states must be evaluated. In general, 
such matrix elements cannot be determined simply, other than by lattice computations. 
However, there are a few special and well known cases where the matrix elements of certain 
gluonic and light quark operators are known by the symmetries of QCD. 

There has been renewed interest in the interactions of the lowest lying hadrons with 
quarkonia. An exciting prospect is measuring the binding of J/ip to nuclei |9|-|TT| 



or even 



measuring the scattering cross section from a single nucleon |[L2| , but theoretical computa- 
tions require knowledge of the couplings between gluons and the J/ip. A large- N c estimate 
of these couplings |||| has been used to predict a binding energy of ~ 10 MeV to infinite 
nuclear matter. Further, the observed anomaly in elastic polarized pp scattering near 5 GeV 



fl3| has been attributed to threshold or resonant effects in cc-N interactions fl4|| . The low- 
momentum interactions of quarkonia are also important in a completely different arena. It is 
possible that the ratio of observed to predicted J / ip production in ultra- relativistic heavy ion 
collisions will provide a "handle" on the formation of a quark-gluon plasma in the collision 
||i"5|j . Whether or not this turns out to be valid, it is important to understand the interaction 



of quarkonia with the debris resulting from such a collision. The existence of co-movers with 
the quarkonia requires that the low-momenta interactions as well as hard interactions 
be understood. 

The natural framework for dealing with low-momentum interactions between the light 
hadrons and quarkonia is chiral perturbation theory. This allows one to relate matrix ele- 
ments for different processes using the underlying symmetries of QCD. Chiral symmetry has 
been used to constrain pionic matrix elements in previous work PJ6Ul7f and further a chiral 
lagrangian in the meson sector was written down at leading order ]TB| , |iT| . This was used to 
determine the leading contribution of the light quark masses to quarkonium level spacings 



U, and to examine the decays ip(2S) -> J/^n + n~, T(2S) T(1S)tt + tt- |J. 

We explore several aspects of the interactions between quarkonia and the low-lying 
hadrons. Firstly we wish to determine if statements can be made concerning the inelas- 
tic scattering of J/ip's by co-moving 7r's as are to be expected in ultra-relativistic heavy ion 
collisions. Secondly, we wish to test the predictions of the large- N c limit of QCD for the coef- 
ficients of the chromo-electric and magnetic operators. These directly impact the predictions 
for the binding of quarkonium to nuclei. Thirdly, we wish to test the usefulness and validity 
of the chirally symmetric lagrange density that matches onto the multipole expansion. To 
this end we explore the radiative decays i>(2S) — > J/ipn+n^j and T(2S) — > Y(lS)n + 7i~ r y. 
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II. THE MULTIPOLE EXPANSION 



The lagrange density describing the interaction between quarkonia and gluons at leading 
order in the gluon momentum and the heavy quarkonium limit is 

C = E (Pt(f)Pv(i) - «) (C^0 E + C^Ob) . (1) 

v QQ 

The subscripts denote chromo-electric and -magnetic interactions respectively, 

O e = -v,v u G» aA G uA 

O b = \g^ a G% - v,v v G» aA Gl A , (2) 

and v a is the four-velocity of the quarkonium. The quarkonia field operators are normalized 
non-relativistically (factors of \/2M removed), where P denotes a pseudo-scalar quarkonium 
and V a vector quarkonium in the same multiplet of heavy-quark spin-symmetry. The % and 
/ that appear as superscripts on the constants Ce,b and as labels on the field operators 
denote the initial and final states of quarkonia, e.g. IS", 2S, etc. By this construction 
we are implicitly assuming that the momentum transfer from the gluons is small, and the 
four- velocity of the quarkonium is conserved during interactions. The mass scale Aqq will 

be of order the inverse radius of the state, and as such we expect the coefficients C^'g to 
be numbers of order unity. In order to progress further in computing observables, these 
coefficients Cg'g need to be determined. As the wavefunctions for color-octet intermediate 
quarkonium states are unknown, one must resort to models or particular limits of QCD for 
theoretical estimates. 

For decays, such as ip(2S) — ► J/ipn + n~ , T(2S) — > T(lS)n + iv~ it was argued that the 
multipole expansion is not applicable as the energy release is the same order as the de- 
nominator in the multipole expansion ||20|| . While it is true that there is no tuneable small 



parameter that can make the multipole expansion arbitrarily accurate for such decays, it is 
possible that numerical factors allow the expansion to be convergent. This might arise from 
the fact that the intermediate states that appear in the matching are in a color octet and will 
not be degenerate with the analogous color singlet states. The scenario appears to work well 
for the strong decays i/j(2S) — > J/ipim and T(2S) — > T(1S)ttt[ as discussed previously and 
in the following sections. However, in order to test this further we will discuss the radiative 
decays ij){2S) — > J/tpTiiT'y and T(2S) — > T(lS)nnj which are calculable at leading order in 
the multipole and chiral expansions. 

The gluonic operators Oe,b can De written in terms of the gluon energy-momentum 
tensor, T^ v 

O b = tyvT + \ GaPAG % . ( 3 ) 



where 

rpfiu _ 

" 4 



= -g^G apA G% - G mA G vA . (4) 
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Higher dimension operators such as those discussed in [[2(J are suppressed by powers of 
the quarkonium radius, or possibly more appropriately the spacing between the color singlet 
and octet states. We will not address the role of the higher dimension operators in the 
multipole expansion. For convenience we will define coefficients Cjjf'g through 

MM 

r (f,i) _ ^E,B /-n 
E ' B ~ A3_ ' 

QQ 

where the units of the C^'b are GeV -3 throughout this work. 

The amplitude for a process, say nJ/ip — > mfj{2S) at leading order in the multipole 
expansion is then given by 

A = {-kiI){2S)\C\-kJ/iI)) 

= (C E + C B ) v,v v (iriP{2S)\TpTrJ/iP) + (C B - C E ) (tt^S^G^G^ttJ/iP) . (6) 



III. MATCHING ONTO THE CHIRAL LAGRANGIAN 

To determine physical scattering amplitudes, matrix elements of T£ u and G a P G^g are 
required between the appropriate light hadron states. For low momentum processes chiral 
perturbation theory is the natural framework to work. The strong dynamics of the octet of 
pseudo-Goldstone bosons is described at leading order by a lagrange density 



2 



f 



where 



d^d^] + ATr[m 9 S + h.c] + ... , (7) 



S = exp(^M) , (8) 



/ 

is the exponential of the meson field (in SU (2) L ® SU {2) R 



"-ex 5 -X,«) • < 9 > 



7T 



and m q is the light quark mass matrix. The dots denote higher dimension operators sup- 
pressed by inverse powers of the chiral symmetry breaking scale, A x ~ lGeV. 



Matrix Elements of G af3A G^ 



It is well known that the divergence of the scale current allows one to determine matrix 
elements of the G a ^ A G A p operator |p|flP, |Il3| , |2~T| at low momentum transfer. The divergence 
of the scale current (the trace of the energy-momentum tensor) for QCD is 



(10) 



4 



where the 7, are the anomalous dimension of the q^i operator and the (3 is the QCD (3- 
function, 



9g 3 



16tt 2 



(11) 



for 3 flavours of light quarks. For our work we will neglect the 7^, as they are small compared 
with unity and work with 



/ 9 i 



(12) 



As the engineering dimensions of the S field of mesons is zero, one finds that the corre- 
sponding quantity for the mesonic part of the chiral lagrangian is 



f 2 

= - J —Tr 



dpWYt] - 4ATr[m q £ + h.c] + ... 



(13) 



where the dots denote terms higher order in the chiral expansion. 

In order to complete the relation between G a/3A G A j3 and operators in terms of the £ 
field, matrix elements of the light quark mass operator are required. Fortunately, this is 
straightforward, 



rriiq^i — > — XTr [m q S + h.c] + 



(14) 



and one finds that 

1 n nvA n A 

2~T G & ^ ' 



P 

1 r 

4 



3\Tr [rriqE + h.c] + 
3 



(15) 



The dots denote higher dimension operators. The matrix element {^{q')\\^-G^ uA G A u \Ti{q)) 
at tree-level is therefore m 2 + p 2 where p = q — q' is the momentum transfer of the it's. This 
result disagrees with |22| but agrees with |J. 

We can compute the matrix element to one-loop level in the chiral lagrangian. The 
higher order counterterms that appear as dots in the previous expression must be retained 
and it is straightforward to show that 

Ag = {<q')\\^G^ A G A Mq)) 
z 9 



ml+p 2 + 



1 



16vr 2 / 2 



~m\p 2 + (3m 2 - 2p 2 )p 2 log I '^f 
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—ml + -m n p — 2p <Px(p /m 



- I0mlp 2 & 2 (p 2 /ml) + wi(fj,)mlp 2 + w 2 {n)p i 
where the functions $1,2 (Z) are defined by 



(16) 
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QJZ) = f 1 dx log (1 - x(l - x)Z - ie) 
Jo 

$ 2 (Z) = / dx x(l -x)log(l -x(l -x)Z-ie) . (17) 
Jo 

The constants are local counterterms arising in the chiral lagrangian whose /i dependence 
exactly cancels the p, dependence of the chiral logarithms. It is natural to choose p = A x , 
the scale of chiral symmetry breaking. We have explicitly retained analytic terms arising 
from the one-loop graph, and have not absorbed them into the counterterms. Since the 
contribution of Wi is always suppressed by factors of m 2 /l67i 2 f 2 compared to the tree- 
level amplitude, we will be unable to fit w\ to data, but we expect that its contribution is 
negligible. It is important to note that the matrix element of G iivA G A v is enhanced by a 
factor of l/a s (p) in its contribution to ipn scattering or nn production. 



B. Matrix Elements of T£ v 



The matrix element of the traceless part of the gluonic component of the energy- 
momentum tensor, T^ u is required in addition to the matrix element of G al3A G A /3 in order to 
determine low-momentum scattering amplitudes. Clearly, we have no symmetry dictating 
the normalization of the matrix elements and we must appeal to data. The operator Tg v is 
traceless and symmetric in its lorentz indices, and as such we match onto an operator in the 
chiral lagrangian at leading order with the same properties, 
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(18) 



where V 2 (p) is a (scale-dependent) constant that must be fit to data. At tree- level in the 
chiral expansion the matrix element of T^ v is 



A? = (n(q')\T g nn( q )) 
= V 2 (jm) ' 



(19) 



and using the variables p = q — q' and k = (q + q')/2 



A 



v 2 (p) 



2 ( /,•"/,•" - { -<rir 



2 \P"P U 



(20) 



This matrix element can be computed simply at one-loop in the chiral expansion. As T^ v 
is traceless in four-dimensions only, we find a contribution proportional to g^ v that yields a 
non-zero trace at one-loop level. This contribution is analytic in the external variables and 
masses and is accompanied by contributions from local counterterms that are unknown. We 
find that 



2 ( Wh v - -g^k 2 



+ 



1 

2 + 



1 



///; 



(2p 2 -m 2 7T ) -log -7 + $ 2(P7<) 
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P^P V ~ -^ V P 2 



6 



16tt 2 / 2 V 6 n 12 ny 6 1 J y 



(21) 



The contributions of and are small and will not be able to be fixed by the available 
data. 

For the processes considered here the scale-dependent constant V^(/i) is not directly 
observable as it always appears multiplied by C^ 1S ^'^ 2S ^ ■ We will use V 2 (A X ) = | 



flqj23|j, however its actual value will not impact this work. Deviations from this value will 
be absorbed into the definition of C^ S ^^ 2S ^ _|_ (ji( ls )'( 2S )\ When a comparison is attempted 
between the extracted values of Cg and (7^ 15 - ) ^ 2 ' s '' ) l with, say, the large- Af c limit of 

QCD, then the precise value of V^(A X ) will have consequence. 



IV. THE STRONG DECAYS ip(2S) J/^vr+vr- AND T(2S) -> T(lS>+vr- 
The matrix element for strong decays ip(2S) — > J '/ 'ifin^ and T(25) — > Y(lS)7r + 7r~ , 

e.g. 

(J/^7t + (q + )n-(q^)\C mt .\ij(2S)) , (22) 

can be obtained from Ac and Ai^ by crossing symmetry. If q + and q_ denote the outgoing 
momentum of the 7r + and 7r~ respectively then we have p = q + + g_ and k = \(q+ — <?-), in 
the expressions for Aq and Atp. 

The data we use are from decays between states in quarkonia to a pion pair. We have 
chosen to analyze the data for T(2S) — > T(lS)7i + n~ p4j and ip(2S) — > J/tpn+n^ |25| and 



not the decays from the T(35) [^6[ as we expect the multipole expansion to be more reliable 
for transitions between low-lying states in quarkonia (for work on strong decays from the 
T(3S') see P7[| ). The data is presented as a normalized differential decay distribution plotted 
against the variable x defined by 

M^s) - M J/4 , - 2m n ' 

for i/j(2S) — > J /iI)'k + 'k~ as shown in fig.l, and an analogous definition for the decay T(2S) — > 
T(lS')7r + 7r _ , as shown in fig. 2. Js^-n is the invariant mass of the tttc system. Note that the 



data we have shown in figs. 1 and 2 is taken directly from fig. 7 of reference |24] and is not 
the original data set. 

The partial-width for each decay modes is taken from the branching fractions and total 
width measurements of the 2S states given in particle data group evaluation |28| 



r(^(2S) -> J/^tt+tt-) = 90 ± 10 keV 
r(T(25) -> T(15)7r + 7r-) = 8.1 ± 1.3 keV . (24) 





s 



Y(2S) J/V|/ 



/X 

✓ v 

FIG. 3. Tree-level diagrams contributing to ip(2S) — ► J/-07r + 7r~. The light-shaded circle de- 
notes insertions of G 2 and T^. 



In what follows we will use only the central values of the above widths and neglect the 
uncertainties. 

At leading order in the multipole expansion, and tree-level in the chiral lagrangian we 
can fit the parameters Cg 1S ^^ 2S ^ and C^ 1S ^^ 2S ^ (from now on we will drop the superscript, 
and simply use Cjj?\ C^ for the charmonium decays and for the T decays) 

to the spectra shown in fig. 1 and fig. 2. The tree-level graphs in the chiral lagrangian 
that contribute to the decay are shown in fig. 3. Using a value for the strong coupling of 
a s (^cc) = 0.55 for charmonium and ot s {r b -^ = 0-33 for bottomonium, we find (V^(A X ) = |) 

- = 7.0 ± 0.2 GeV- 3 , Cf + = 14 ± 1 GeV~ 3 
CP - C { P = 1.3 ± 0.1 GeV" 3 , Cg° + C { P = 3.9 ± 0.7 GeV" 3 . (25) 

In fitting the data we have weighted each data point equally (with an error of ±0.1) and 
not used the experimental uncertainty indicated by fig.l and fig. 2. We find a much tighter 
constraint on the difference between Cb and Ce (the coefficient of the G 2 operator) than on 
the sum (the coefficient of the T^ v operator). These extractions then lead to 

eg = 10.5 ± 0.6 GeV" 3 , C^ = 3.5 ± 0.5 GeV" 3 

C { P = 2.6 ± 0.4 GeV" 3 , C [ p = 1.3 ± 0.3 GeV" 3 . (26) 

We note that in fitting the data there is a global sign ambiguity in the coefficients Ce,b- 
We have chosen Ce > for our analysis. The uncertainties are from the fit only and do 
not include the uncertainty associated with the widths. We see several interesting trends in 
these extractions. Firstly, the chromo- magnetic coefficient Cb is significantly smaller than 
the chromo-electric coefficient Ce for both systems. This is expected on grounds that the 
magnetic interaction should be suppressed by a power of the heavy quark three- velocity in 
the quarkonium. Secondly, both coefficients in the T system are substantially smaller than 
the corresponding quantity in the charmonium system. Again, this is to be expected on 
the grounds that the strong interactions will decouple in the limit of vanishing radius of the 
system. Our results at this order do not agree with results found in [|TJ]] where coefficients 
found in that work were approximately the same for both systems. 

Fitting at one-loop is complicated by the fact that there are seven counterterms W\ —Wj. 
The counterterms w\, w 3 , w 5 , Wq have only a small impact upon the decays under consider- 
ation as they are suppressed by powers of the 7r mass compared to tree-level terms. The 
counterterm wj occurs in the trace of T^ v at loop- level and as we have neglected the trace of 
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FIG. 4. Tree-level and loop-level coefficients C E and C B for ^(25) - J/ijj and T(25) - T(15) 
systems. The 1 — a regions arise from a Monte-Carlo search in parameter space. We have used 



V 2 (A 



2 ' ^ s 



0.55 and a s {r b -r) = 0.33. The loop-level fit is for = 0. 



T^ v arising at loop- level (it is analytic in momenta and masses) we have set w-j = 0. Naive 
dimensional analysis suggests that all the counterterms Wi are of order unity, and hence the 
effects of setting each Wi = should be small. However, we find this not to be the case. It 
is clear from figs. 1 and 2 that the decays are well described at tree-level and additional 
momentum dependence arising at loop level will be small. As the amplitudes are dominated 
by the G 2 operator, we suspect and find that a non-zero value for the w 2 counterterm is 
required. With so many counterterms we can do little more than give examples. 

As an explicit example, we set = and have w 2 as the only counterterm to be fit to 
data. If the pion energy distributions became available then we would be able to fit both 
102 and 104, but at present this is not possible. The counterterms w\-w-j are the same for all 
quarkonia systems, and we therefore simultaneously fit Q , Ce\ and w 2 to the 
data sets for ip(2S) — ► J/ipn+'K^ and T(2S) — > T(lS)n + 7T~ . The loop-level graphs in the 
chiral lagrangian that contribute to the decay are shown in fig. 5. We find that 

= 7.8 ± 0.4 GeV" 3 , = 1.3 ± 0.3 GeV~ 3 

= 2.0 ± 0.2 GeV~ 3 , Cjj° = 0.75 ± 0.14 GeV" 3 
w 2 (A x ) = -6.9 ±0.4 , (27) 

where the errors on each quantity are correlated. It is comforting to see that the value 
of coefficients extracted at one-loop level are only slightly reduced from their tree-level 
extractions. One might fear that the relatively large momentum transfer involved in these 
processes ~ 590 MeV may make the chiral expansion invalid, but it would appear that this 
is in fact not the case. However, there is a hint of a problem by the large value of w 2 . 
The value of w 2 is such that it almost cancels the 0(p 4 ) contribution from the loop. It is 
interesting to note that for this fit we find that Cb is much smaller than Ce in the -^-system. 
We are able to fit both the ip and T spectra with single counterterm that is independent of 
the quarkonium system. This suggests that higher order terms in the multipole expansion 
are small. 

Fitting both w 2 and leads to a slightly better fit but the uncertainties associated 
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FIG. 5. One- loop diagrams contributing to ip(2S) — > J/ifjTT + 7r . The light-shaded circle denotes 
insertions of G 2 and Tg v . The dark circle denotes a strong interaction in the chiral lagrangian. 

with the parameters are large. The central values for the Ce,b are very different from the 
tree-level extractions and one might worry that the expansion is breaking down, but the 
large uncertanties preclude any firm conclusions. 

We note that others have looked at the role of final state interactions in these decays 
ZS|], but not using the framework of chiral perturbation theory. 



V. COMPARISON WITH THE LARGE- A^c* LIMIT OF QCD 

One is presently unable to match directly from QCD onto the lagrange density involving 
quarkonia and gluons. Previous estimates of the interaction between quarkonia and light 
hadrons have been based largely on estimates of the various Ce and Cb in the large- N c 
limit of QCD |4]|J. In the large- N c limit the repulsion in the colour octet channel vanishes 
as 1/N C while the attraction in the singlet channel remains. Consequently, the colour octet 
intermediate state in the dominant diagrams in matching onto Ce and Cb can be described 
by plane waves. Further, the external color singlet states are coulombic. Matrix elements 
of a double insertion of the chromo-electric dipole operator are simply computed in terms 
of the overlap between coulomb bound-state wavefunctions and plane waves with a plane 
wave propogator. For zero-momentum transfer amplitudes this gives a sensible estimate for 
the actual matrix element, as there is a mass-gap between the intermediate state and the 
bound-state. However, for decay processes such as i>(2S) — > J '/ipim, the estimate is assured 
to be unreliable as the energy release is comparable to the gap between the bound state and 
the intermediate state(s). 

At leading order there is no contribution to Cb while there is a calculable contribu- 
tion to Ce- For dipole-dipole interactions between the same initial and final states it is 
straightforward to show that 

r [(nS),(nS)] _ 1 3 j 

^e — 2 o n 
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(i is summed over) where ao is the bohr radius of the quarkonium ground state Oq ~ 
(0.7 GeV) -1 , aj ~ (1.5 GeV) -1 f30fl . The (p n are coulomb wavefunctions for the n-th S- 



state, e is the binding energy of the n — th state and eo is the binding energy of the ground 
state. This calculation leads to 

_ n s U7T r [(is),(is)] _ n 
— o~27~ ' B ~ 

C K2SU2 S)] = a 3l0047T C i(2S),(2S )]=Q (2Q) 

— I 

which gives 

C [{isum* „ 4.7 (GeV )-3 , C fsws)W „ 341 (GeV) -3 

C K15),(15)]T ^ _ 4g (GeV) _ 3 ) C K25) )( 25)]T ^ 35 (3()) 

There is a problem in estimating the C^ 2 J' ) ' ( ' 15 ' ) ' from this construction from the fact that 
the energy release in the transition is the same magnitude as the off-shellness of the inter- 
mediate states in the free colour-octet propogator. In this limit a short- distance expansion 
makes little sense |2(| and the operator-product expansion we have employed is inappro- 
priate. However, it is more than likely that the colour-octet intermediate state has a mass 
splitting to the colour- singlet states, although the splitting does not become large in any 
limit of QCD. The splitting would allow one to perform the OPE as we have done and for 
the operators Oe,b to be the leading terms in a systematic expansion. As such, we cannot 
directly use the construction of the large- N c limit of QCD to make a reliable estimate of 
the C^ 2 ^'^ 1S ^. We can compare numerically the values of C^ 2 J^ 15 ^ we have obtained in 
previous sections and look to see if at least they are consistent. If we naively assume 



W(2S),(1S)] _ C [(1S),(1S)] C [(2S),(2S)] 

to provide an extremely rough estimate of what we might expect for the off-diagonal element, 
then we find that C [ ^ 2SU1SM ~ 40(GeV)~ 3 and c^ 2S) ^ ls)]r ~ 4(GeV)~ 3 in the large-iV c limit. 
These estimates are within factors of ~ 4 and ~ 2 of our extraction. From this we very 
cautiously hope that the large- N c estimates of the coefficients C^ n ^'^ nS ^ are reasonably close 
to their true values and that the estimates for the binding energy of quarkonia in nuclear 
matter [|H|] are not unreasonable. 



VI. APPLICATIONS 

A. Inelastic irJ/ip -> tti(j(2S) and ttT(IS) -> vrT(25') Scattering 

With the possibility of " suppression" being a tool for detecting the formation of a 
quark-gluon plasma in ultra-relativistic heavy ion collisions it is important to understand 
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FIG. 6. The total cross-section in millibarns for nJ/tp — ► inp(2S) as a function of the tt lab 
energy in GeV. The curves are evaluated at the value of parameters corresponding to the best fit 
to ip(2S) — > J/ipTr + ir~ at tree-level and loop-level (with W2 = —6.9 and the remaining Wi = 0). 

inelastic J/ip scattering in a gas of hadrons. Naively, one expects 7r's to comprise a non- 
negligible fraction of the number density in a thermal hadronic gas, and as such, 7rJ/ijj 
inelastic scattering needs to be understood as a background to other mechanisms that remove 
the J/^'s produced in any given collision. 

For wJ/if) -> 7n/>(2S) and nT(lS) -> ttT(25) near threshold (£[ ab ~ 810 MeV) we can 
use directly the expression for Aq and A? discussed previously. It is likely that there will 
be significant corrections to the amplitudes for irJ/tp — ► mp(2S) and 7rT(15) — > ttT(2S) 
determined from the crossed channels T(2S) — > T(1S)tt + 7t~ and ip(2S) — > J/ip7T + 7r~ simply 
due to the large energy transfer involved and a break-down of the chiral expansion in this 
kinematic regime. We expect that the calculation will fail close to threshold and the differ- 
ence between the cross-section as computed at tree-level and at one-loop level will give some 
feeling for the range of validity of the chiral expansion. The curves in fig. 6 indicate that the 
loop-level and tree-level estimates are not so different even near where the expansion must 
fail. This should be taken with extreme caution. However, we see that the cross section is 
of the order of 10 -2 mb in the region where the computation is most reliable. 

B. Radiative Decays V(2£") -► J/^tt+tt^ and T(25) -» T(lS , )vr+7r-7 

The radiative decays ip(2S) — ► J/tpn+Ti^j and T(25) — > T(lS')7r + 7r _ 7 are dominated 
by electromagnetic interactions of the 7r's that are calculable with the chiral lagrangian out- 
lined in the previous sections, after gauging with respect to the electromagnetic interaction. 
Operators involving the electromagnetic field-strength tensor F^ u arising in the matching 
between QCD and the effective theory of quarkonium will be suppressed by two powers of the 
quarkonium radius. There will be higher dimension operators in the chiral lagrangian arising 
from matching the theory of quarkonia and gluons onto the chiral lagrangian and quarkonia. 
Such operators are suppressed by additional powers of the chiral symmetry breaking scale. 
The leading contributions to ip{2S) — > J /ipn + ix~^ and T(2S) — > T(lS')7r + 7r~7 occur at tree- 
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FIG. 7. The tree-level diagrams contributing to ip{2S) — ► J/ipTr + ir 7. The light shaded circle 

9 



denotes insertions of the G 2 and T£ v operators. 



level and hence we do not need to discuss the role of local counterterms at leading order. 
We note that the decays ip(2S) — > J/ip7T 7T 0r y and T(2S) — > T(15 , )7r°7r°7 are forbidden by 
charge conjugation. 

One might naively expect a significant contribution to the ttti'j decays from pole graphs 
with x c o(lP), Xd(lP) and x c2 (lP) as intermediate states. However, G-parity forbids the 
Xcn(lP) — > J/'ipim transitions and therefore, such contributions are not present. 

It is straightforward to determine the tree-level matrix element for the radiative decays. 
There are contributions from both the G 2 operator and the T£ v operator as shown in fig. 7. 
The G 2 operator yields an amplitude of 



p + ■ e p- ■ e 
k ■ p + k ■ p_ 



(32) 



where k is the photon momentum, is the 7r + ' momentum resp., and e is the photon 
polarization tensor. The contribution from T£ v is found to be 



( v ■ kp + ■ e \ 

V • p_ V ■ £ I — V ■ p + 

\ p+ ■ k J 




%R T = -2ieV 2 (fx) 

+ (v ■ p + v ■ p_ - ^p + • p_ - ifc • (p+ + p-) ) ( 7^-3— - 77- — ) • (33) 

The full amplitude at this order is given by 

%A 1 = (C E + C B ) \in T ] + (C B - C E ) \iK G ] . (34) 

The distribution dT/dY resulting from %A^ as a function of the normalized photon energy 
Y = E-f/E™*- (P™ ax - = 309 MeV) is shown for different values of the coefficients C E and 
Cb in fig. 8. In fig.9 we have shown the same differential spectrum weighted with the square 
of the normalized photon energy to highlight the higher energy part of the distribution. 
Experimentally verifying the differential distribution will be a good test of the chiral and 
multipole expansions. 
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0.2 0.4 0.6 0.8 y 10 

FIG. 8. Differential decay spectrum as a function of normalized photon energy in 
ip(2S) — ► J / 'ipir + tt~~ 7 . Shown are the distributions obtained for different values of Ce and Cb- 
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Y ■ 

1 mm 


Branching Fraction 


0.2 


1.0 x 10~ 3 


0.4 


2.2 x 10" 4 


0.6 


3.5 x 10" 5 


0.8 


1.8 x 10~ 6 



TABLE I. Branching fractions for ip(2S) — ► J/ifjTT + ir 7 for different cuts on the normalized 
photon energy, Y min , for a s {r cZ ) = 0.55, C E = 10.8 GeV~ 3 and C B = 3.8 GeV~ 3 . 

The integrated branching fraction is infrared divergent as we have not included one-loop 
vertex and wavefunction graphs. However, it is useful to have the branching fraction for 
events with a E 1 greater than some lower energy E™ m . In table 1 we have given the branching 
fraction for different values of Y min for the tree-level extracted values of Ce = 10.8 GeV~ 3 
and C B = 3.8 GeV~ 3 . 



C. Electromagnetic Polarizability of Quarkonia 

While it is more than likely that the electromagnetic polarizability of quarkonia will never 
be measured, there is an interesting theoretical point to be discussed. The leading contri- 
bution to both the electric and magnetic polarizabilities, defined by an effective lagrange 
density, 

£ eS - = E (Pl(f)Pv® - t£,(/)K"(0) [~Ma + P)v<rV p F° s F s p + Tr(3F° s F aS \ , (35) 

V 

comes from tree-level pole graphs, where two photons are coupled to the quarkonium via 
colour-singlet intermediate states. The constants a and (3 are the electric and magnetic 
susceptibilities. Such contributions are finite in the chiral limit and are naively estimated 
to be 

a ~ 3 — r 3 - (36) 

in the same way that Ce,b are estimated in the large- N c limit. However, there is a contri- 
bution to a and (5 from Ce and Cb via n loops coupled to two photons, as shown in fig. 
10. This is numerically subleading to the contribution from the pole graphs, except in the 
extreme chiral limit. In this limit the loop graphs are logarithmically infrared divergent, 

a^ = ^ = ^VM(C E + C B )^lo g (^) , (37) 

and only Tg v contributes. Using the large-iV c estimates of Ce = 4.7 GeV -3 and Cb = 
GeV~ 3 along with the actual value of the ir mass and V^/i) = |, one finds a loop contri- 
bution of a loop = (3 loop ~ lO" 6 fm 3 . 
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FIG. 10. Diagrams contributing to the electromagnetic polarizabilities of the J/tp. The shaded 
circle denotes insertions of G 2 and Tg V ■ 



D. Decays from the T(35) 

Until this point we have been considering decays from the 2S states in the ip an d T 
systems. However, there is data available for the transitions T(35) — > T(1S)7T7t as shown 
in fig. 11, and T(3S) — ► T(2S)mr |26|]. (Note that we have taken the data from fig. 11 



of [ 26 1 and it is not the original data set.) The spectrum for T(35) — > T(2S)tctt is similar 
to that of the decays from the 2S states and we will not discuss it further. There has 
been significant discussion of T(35) — > T(lS)im as the differential spectrum does not have 



the same shape as that of transitions from the 2S states p9|j31| , |32| . It has been suggested 
that a resonance in the T(1S)tt channel P9"| , |3"T| is responsible for the observed distribution. 
Another suggestion was that the multipole expansion is subdominant to hadronic effects 



involving B-mesons |B2|. It is important to note that the invariant mass of the irn system 



can extend up to 0.9 GeV. The soft pion analysis is expected to fail for the higher mass 
pairs. As such, one might expect to see modifications to the differential spectrum before 
the end-point. Further, it was assumed that Cb vanished in the previous analyses, allowing 
only for a spectrum peaked near the high end of the spectrum. It is clear from fig. 11 that 
the best tree-level fit to the entire range of x is poor, with the apparent minimum at x ~ 0.6 
not reproduced. If we restrict ourselves to the low end of the spectrum where the chiral 
lagrangian is expected to be applicable, then that portion of the spectrum can be described 
with the multipole expansion. T^ u dominates the decay amplitude with Ce ~ Cb ~ 0.3 
yielding the best fit for m nn < 0.65GeV. Clearly, only a portion of the spectrum can be 
described at leading order in the multipole and chiral expansions. Whether the deviations 
from the leading result arises from higher order terms in the chiral or multipole expansions 
or from a given resonance cannot be determined. However, in order to reproduce the low- 
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mass region of the spectrum, a large value for the chro mo-magnetic coupling Cb is required, 
in contrast to transitions from the 25* states. 



VII. DISCUSSION 



We have explored the strong interactions and radiative decays of quarkonia in order to 
better understand the multipole and chiral expansions that might be appropriate for such 
processes. The unproven convergence of the multipole expansion is assumed throughout 
this work and we suggest that numerical factors could lead to convergence despite there 
being no rigorous limit of QCD in which this is true. Current discussions of how quarkonia 
interact with the light hadrons away from the chiral limit have been addressed and the 
strong decays ij)(2S) — > J/tpir + ir~, T(2S) — > T(lS)7r + 7i~ were focused on at tree-level 
and loop-level in the chiral expansion. The leading couplings in the multipole expansion 
are found to be significantly smaller in the T system than in the ip system, consistent 
with naive scaling arguments. Further, counterterms appearing at loop level in the chiral 
expansion are found to be approximately the same for the ip and T systems, suggesting 
that the multipole expansion is converging. In addition, we discussed the large-iVc* limit 
in which estimates of the coefficients Ce,b f° r various quarkonium states have been made. 
While there is nothing concrete that can be concluded, it would appear that the large- Nq 
limit is not in serious conflict with data. Crossing symmetry has been used to estimate 
the size of the inelastic nJ/ip — > nip(2S) scattering cross section that might be relevant 
for understanding ip suppression in ultra-relativistic heavy ion collisions. In the region 
where the theory is expected to be most reliable we find a small cross section, ~ 10~ 2 mb. 
Power counting in the quarkonium radius indicates that the radiative decays ij){2S) — > 
J/'07r + 7r~7, T(2S) — > T(lS)ir + TT~ / -f are dominated by tree graphs computable in chiral 
perturbation theory. The photon energy spectrum is presented and branching fractions 
~ 10~ 4 are found. These branching fraction depends sensitively upon the values of Ce,b 
used and thus provides a consistency test of the calculational framework. Finally, we made 
brief remarks on transitions from the T(3S) and on the electromagnetic polarizability of the 
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